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I. INTRODUCTION 

It is well-known that fourth derivative terms are nec¬ 
essary in the gravitational action, to render the theory 
consistent at the semiclassical level (see mm for an intro¬ 
duction to the subject and also [3] for a recent pedagog¬ 
ical review). In Quantum Gravity (QG) the same fourth 
derivative terms provide renormalizability of the theory 
gj. At the same time, adding fourth derivative terms to 
the standard Einstein-Hilbert action has an undesirable 
effect, since the spectrum of the theory in the spin -2 sec¬ 
tor gains a massive unphysical ghost in addition to the 
usual healthy graviton. An attempt to remove these un¬ 
physical states from the initial quantum state leads to 
a non-unitary quantum theory. In the semiclassical the¬ 
ory, when gravity is a classical background for quantum 
matter fields, there is no problem with the unitarity of 
the gravitational S'-matrix, hence the condition of con¬ 
sistency can be reduced to the requirement of stability 
of the physically relevant classical solutions with respect 
to small metric perturbations. In the present work we 
will discuss a new aspect of this problem, related to the 
difference between the massive unphysical ghost and the 
tachyon. 

Let us start by fixing the notations. The action of 
gravity corresponding to the renormalizable semiclassical 
theory (for introduction one can see, e.g., [IH3], further 
references therein and also recent more formal work 0) 
consists of the following terms: 

Svac = Seh + Shd , (1) 
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where Seh is the Einstein-Hilbert action with a cosmo¬ 
logical constant term, 


Seh = ~\Lg / dAx ^g{R + m 


( 2 ) 


and Shd includes a minimal necessary set of higher 
derivative terms, 

Shd — J d^X\/ —g{a,\C 2 + 02^4 + a 3 lUI? + 04 !?“}.(3) 


In this expression 


£4 — R^ap — 4i ? 2 


a/3 


R 2 


( 4 ) 


is the Gauss-Bonnet term (Euler density in d = 4) and 
C 2 is the square of the Weyl tensor. Furthermore, we 
use units corresponding to c = h = 1 and denote with 
M p = 1/VG the Planck mass. The metric has signature 
(+, —, —, —). For the sake of simplicity we assume that 
the absolute value of ai is of order one. 

A general observation about the form of the action 
and especially higher derivative terms ([3]) is in order. 
One can deal with these terms in two different (albeit 
physically equivalent) ways. First, it is possible to con¬ 
sider the classical action of gravity being free of higher 
derivative terms. Then one can observe that (different 
from the flat-space QED, for instance) there are higher 
derivative divergences, corresponding to the running of 
these terms and hence the non-local structures such as 


we shall discuss later on in Eq. (22). This means that 
even if the terms ([3]) are not introduced at the classical 
level, the corresponding quantum contributions emerge 
anyway. Second, one can follow the standard QFT ap¬ 
proach and introduce into the classical action all those 
terms which will emerge with divergent coefficients in the 
loop corrections. The main advantage of this procedure 
is technical simplicity, because one can always deal with a 
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usual renormalization procedure, consider the renormal¬ 
ization group running in a useful conventional way [2] etc. 
The last observation is that the list of the terms in ©> 
and © follows from the fact that only these structures 
emerge in the semiclassical approach at all loop orders. 
This fact is known starting from the classical work of 
Utiyama and DeWitt [B], and can be easily seen from the 
covariant renormalizability and power counting for the 
semiclassical theory 0. 

The spectrum of the theory ([!]) linearized around a 
Minkowski background contains in the spin-2 sector the 
massless graviton and an additional spin-2 particle with 
mass of order M p . If or is negative this particle has 
positive mass square, negative kinetic energy and hence 
is a ghost. For a positive ai the mass square is negative 
and the particle is a ghost and a tachyon at the same 
time [7]. We postpone the detailed description of these 
two types of particles to the next section and now only 
make some general observations. 

According to the recent discussion in 0, at least on a 
cosmological FRW background massive ghost excitations 
are not generated far below the Planck threshold of fre¬ 
quencies. The same result was found much earlier PM 
for de Sitter background. Qualitatively similar consider¬ 
ations can be found in m for generic models of mas¬ 
sive gravity, where one can also meet massive ghosts or 
tachyons. The result of [5] means that a massive unphys¬ 
ical ghost can be destructive for the theory, but only if it 
appears as a real physical particle. On the other side, the 
presence of a ghost in the mass spectrum of the theory 
may not affect stability if such a ghost is just a virtual 
particle in the vacuum state of the theory. One can note 
that the problem of massive ghost is a tree-level prob¬ 
lem of quantum theory. It is equivalent to the absence 
of stability in the classical theory. In practical terms the 
conclusion of |8j means that one can choose the initial 
condition in such a way that the system does not de¬ 
velop instability for sufficiently long time. According to 
[8] this choice of initial conditions is possible for ai < 0. 

On the contrary, for the case a± > 0 the extra mas¬ 
sive particle is a ghost and tachyon at the same time and 
satisfies an anti-oscillator equation with exponential so¬ 
lutions. Then instabilities are unavoidable, even without 
the presence of an external force. Moreover, the intensity 
of the instability may be even enhanced if the mass of the 
tachyon increases. There is no possible choice of initial 
conditions which can help to control the instabilities in 
this case. 

One can naturally ask whether we can define ai to be 
zero and thus avoid to have both ghost and tachyon states 
in the spectrum. The answer to this question is negative, 
because at the quantum level the parameter ai runs due 
to the change of the energy scale. In this work we will 
assume that gravity is classical and we take into account 
only quantum effects of matter fields. In the last section 
we will also comment on the possible role of quantum 
gravity within some of the existing approaches. 

Let us stress that the running of ai is a pretty well- 


explored theoretical phenomenon. In particular, the full 
expressions for the “physical” /3-functions in the momen¬ 
tum subtraction scheme of renormalization were calcu¬ 
lated in H2 and [T^|. At low energies one can observe 
quantum decoupling, which is a close analogous of the 
Appelquist and Carazzone theorem in QED [16j. This 
means that in the present-day universe, when the typ¬ 
ical energy scale of the background is of the order of 
Hq ~ 10 -42 GeV , all massive fields are too massive to 
affect the running of a\ in the IR. Hence, in the late uni¬ 
verse it is sufficient to take into account quantum effects 
produced by the unique massless particle, which is the 
photon. Starting from Sect. Ill it will be done in the 
framework of anomaly-induced effective action, which is 
the most appropriate formalism for classically conformal 
fields. 

It is well-known that the renormalization group in 
gravity meets a hard problem in identifying the energy 
scale /x with some physical quantity related to gravity 
(see, e.g., the discussion of this issue in na and also 
the scale-setting procedure suggested in mm)- At the 
same time, all these difficulties concern only quantum ef¬ 
fects of massive fields, while the case of photon is differ¬ 
ent. The point is that electromagnetic fields are massless 
and possess local conformal invariance. Thanks to this 
symmetry one can integrate the conformal anomaly EDI 
[22] and arrive at the reliable form of the semiclassical 
corrections to the classical action. The anomaly-induced 
effective action [23] can be regarded as a local version of 
renormalization group [23 j. 

The purpose of the present work is to use the anomaly- 
induced action to explore the running of or in the far IR 
due to loop effects of virtual photons on a cosmological 
background. This running predicts that at some point 
in the distant future the parameter ai changes sign and 
becomes positive. At this instant the tachyonic mode 
with mass parameter unbounded from above m 2 ~ 1/ai 
will manifest an explosive growth of gravitational waves 
within the full range of frequencies starting from zero 
to the Planck scale. This means that in a very short 
time interval metric perturbations break down the linear 
regime. 

The paper is organized as follows. Sect. [TT] mainly 
serves pedagogical purposes and describes the classifica¬ 
tion of ghosts and tachyons in free theories with actions of 
both second- and fourth-order in derivatives. We also dis¬ 
cuss of the difference between ghosts and tachyons in the 
effective field theory approach. In Sect. m we summarize 
the main results concerning the anomaly-induced effec¬ 
tive action of gravity and explore the effective equations 
of motion on cosmological backgrounds for radiation-, 
matter- and cosmological constant - dominated solutions. 
Sect. |IV|includes the final results for the physical running 
of ai in the same three cases. In Sect. [V]we specialize 
our study to the late de Sitter phase of the evolution of 
the universe and estimate the time period of time for the 
universe to accelerate until the instant when a\ changes 
sign from negative to positive. Finally, in Sect. |VI| we 
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critically discuss the list of assumptions which have been 
introduced in our analysis and draw our conclusions. 

II. GHOSTS AND TACHYONS 

In this section we briefly describe the classification 
of free fields into normal “healthy” ones, ghosts and 
tachyons. For our purpose it is sufficient to consider flat 
space-time. 

A. Second-order theories 

For the second-order theory the general action of a free 
field h(x) = h(t, r) is 

S(h) = J d A x {ri^d^hd^h — S 2 'm 2 h 2 } 

= j d 4 x {ti 2 — (\7h) 2 — S2m 2 h 2 } . ( 5 ) 

Both si and S 2 are sign factors which take values ±1 for 
different types of fields. In what follows we consider all 
four combinations of these signs. 

It proves useful to perform the Fourier transform in 
the space variables, 

h(t,r) = J d 3 k e lk r h(t, k), (6) 

and consider the dynamics of each component h = 
h(t, k) separately. It is easy to see that this dynamics 
is defined by the action 

S'k(Zi) = J dt {A 2 — k 2 h — S2m 2 h 2 } 

= y/ dt{h 2 -m 2 k h 2 }, (7) 

where 

k 2 = k ■ k , and m k = s 2 m 2 + k 2 . (8) 

The properties of the field are defined by the sign of Si 
and s 2 - The possible options can be classified as follows. 

(i) Normal healthy field corresponds to si = s 2 = 1- 
The kinetic energy of the field is positive and the minimal 
action can be achieved for a static configuration. Also, 
the equation of motion is of the oscillatory type, 

h + m k h = 0 , (9) 

with the usual periodic solution. 

(ii) Tachyon has Si = 1 and s 2 = —1. The classical 
dynamics of tachyons is described in the literature, e.g., 
[251 126 ] . hence we present just a very brief review. 

For relatively small momenta m 2 < 0 in Eq. (jsj) and 
the equation of motion is 

h — uj 2 h = 0, u) 2 = \m k \ . (10) 


Then we have an anti-oscillatory equation with exponen¬ 
tial solutions, 

h = h ie ut + h 2 e~ ult , uj 2 = \m 2 k \ . (11) 

If the particle moves faster than light the solution is of the 
oscillatory kind, indicating that such a motion is “natu¬ 
ral” for this kind of particle [25, ' ! 26j. 

An additional observation concerning tachyons is in or¬ 
der. For a field interacting with an external gravitational 
background, it is not impossible to have a situation where 
the same wave changes from being a normal healthy state 
to a tachyonic one in different parts of the space-time 
manifold. In principle, this situation may produce very 
strong effects at both quantum m and classical [28] lev¬ 
els. In the present paper we are considering a qualita¬ 
tively similar situation, where the massive component of 
the gravitational wave changes its fundamental proper¬ 
ties due to the change of the energy scale or due to the 
time evolution in de Sitter vacuum. 

(in) Massive ghost has Si = —1 and s 2 = 1- It is not 
a tachyon, because m 2 > 0. In this case the kinetic en¬ 
ergy of the field is negative, but one can postulate zero 
variation of the action and arrive at the normal oscilla¬ 
tory equation A particle with negative kinetic energy 
has the tendency to achieve a maximal speed, but a free 
particle can not accelerate, for this would violate energy 
conservation. Hence a free ghost does not produce any 
harm to the environment, being isolated from it. How¬ 
ever, if we admit an interaction with healthy fields, the 
tendency of a ghost is to accelerate and transmit posi¬ 
tive energy difference to these healthy fields |29j, e.g., in 
the form of quantum emission of the corresponding par¬ 
ticles. There are many detailed discussions of the fun¬ 
damental problems at both classical and quantum levels 
which arise in the theories with ghosts (see, e.g., 30] and 
EH). In general, the problem of ghosts and the con¬ 
sequent conflict between renormalizability and unitarity 
is one of the most important for quantum gravity, and 
hence attracted a lot of attention. The list of references 
on the possible approaches to avoid the problem of higher 
derivative ghosts can be found in Ref. [8]. In this work 
we also suggested a new approach to dealing with ghosts, 
which is based on the effective field theory ideas, but is 
technically classical and very simple. There are strong in¬ 
dications that the creation of ghost states from vacuum 
is strongly suppressed for a relatively weak cosmological 
background of the late universe. The consideration pre¬ 
sented below is based on the same idea, and assumes the 
same approach. 

(iv) Tachyonic ghost has Si = s 2 = —1. For relatively 
small k 2 we have m k < 0. The kinetic energy is negative 
and the derivation of the equations of motion requires an 
additional definition similar to the non-tachyonic ghost 
case. After that, one can notice that the equation of 
motion is of the anti-oscillatory type, eq. ( [l0| ) and the 
solutions are exponential 
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B. Fourth-order gravity at the linearized level 


In the fourth-order gravity (|T|) the equations for the 
metric perturbations in flat space can be easily obtained 
from the more general ones on cosmological backgrounds 
m (see also Eq. ( [52]) in Sect. 0, 

h + 2k r h + k i h - 1 - (h + k 2 h)=0. ( 12 ) 

327rGai v ' 

It proves useful to introduce a new notation 

1 2 
= - s 2 m , 


327rGai 


(13) 


where S 2 = —signal and to 2 > 0. Then one can recast 
Eq. (12) into the form 


(14) 


where = k 2 + S 2 m 2 . The solutions of the last equa¬ 
tion can be different, depending on the sign of ai and 
hence of S 2 . The general formula for the frequencies is 


(d 2 , 2 \ 

( d 2 2 \ 

\W +k ) 

\W- +rrik ) 


X / 2 

UJ 12 ~ ±*(k 2 ) and ui 3 ^ ~ ±( — to 2 ) 


2\-i/2 


(15) 


For a negative ai there are only imaginary frequencies 
and hence oscillator-type solutions. On the contrary, for 
a positive ai we have S 2 = — 1 and the roots 1^4 are 
real, since in this case — to); > 0 for sufficiently small k 2 . 
Indeed, the first couple of roots corresponds to the mass¬ 
less graviton, and the second couple to an extra massive 
particle. According to our classification, this particle is a 
ghost for ai < 0 and, simultaneously, ghost and tachyon 
for ai > 0 (see [7] for a detailed discussion). 


C. Ghost vs tachyon 


there are strong indications that the gravitational the¬ 
ory with massive ghosts can be free of instabilities at low 
energies, since for £ <C Mp the ghost may be a virtual 
particle and effectively it is not created from vacuum (Sj 
(including Erratum of this paper). The reason is that a 
weak gravitational background is not providing sufficient 
density of energy to generate ghost as a real particle ex¬ 
citation. Only for the typical frequencies of the Planck 
order of magnitude ghost becomes destructive, while at 
much lower energies one does not need to worry about 
its presence in the spectrum. By the end of the day the 
situation is very close to the one in effective field theories 

E2- 

On the contrary, no low-energy protection can be ex¬ 
pected in the theory with tachyons, because they produce 
instabilities independently on their interaction to nor¬ 
mal particles or on the intensity of the background. In 
other words, for tachyons the exponential behavior 
occurs at all frequencies, and not only above the Planck 
threshold [5]. Therefore, the difference between ghosts 
and tachyons is expected to be critical for the low-energy 
regimes. 


III. ANOMALY-INDUCED EFFECTIVE 

ACTION IN THE COSMOLOGICAL SETTING 

In what follows we will use the formalism of anomaly- 
induced effective action of vacuum, so let us start by 
describing an application of this approach to the late- 
epoch cosmology. 

For a general theory including N s massless conformal 
scalars, Nf massless Dirac fermions and N v massless vec¬ 
tors, the anomalous trace of the energy-momentum ten¬ 
sor is [T] 

<T£> = - (/3iG 2 + P 2 E 4 + foOR ), (16) 


The main difference between ghosts and tachyons is 
that a ghost may cause instabilities only when it couples 
to some healthy fields or to the background, while with 
tachyons there is no such a protection. 

The situation with ghosts can be kept under control 
in the effective field theory framework, as it was recently 
discussed in [32]. In an effective field theory there may 
be an apparent ghost due to the low-energy expansion 
in the powers of E/M , where £ is the energy and M 
is cut-off. At the same time, there could be no ghost 
in the underlying fundamental quantum theory. In this 
situation one can safely use an effective description for 
energies £ <C M and do not pay much attention to the 
presence of ghosts. The situation in gravity is similar, 
but with an important difference that the ghost-free UV 
completion of the renormalizable theory is not known. 
Let us note that string theory can not be regarded as 
such an UV completion, since the -R 2 „-like terms which 
are the source of ghosts are removed in string theory by 
a special transformation of the background metric [33] 
and not by a low-energy expansion. On the other side, 


( 47 0 2 & = t 4 , N s + ^ N f + ^N v , 


10 

( 4 tt ) 2 fa = - — N a - — N f - — N v , 
v ’ M 360 360 ; 180 

( 47 r ) 2 & = — N s + — N f - — N v . 
v ; M3 iso 30 2 10 


(17) 


It was already mentioned in the Introduction that in 
the late universe the loops of all massive fields decou¬ 
ple from gravity. The reason is that the typical energy 
scale of the universe is usually defined by the Hubble pa¬ 
rameter, which has the present-day order of magnitude 
Hq oc 10 -42 GeV. This is a very small value if compared 
even to the lightest known particles, e.g., the standard 
estimate for neutrino is to„ oc 10 -12 GeV . Therefore, 
according to the existing results on the gravitational de¬ 
coupling nails] in the present-day (and certainly later) 
universe one need to take into account only the contri¬ 
bution of photons, which are massless. Hence one has to 
set N s = Nf = 0 and N v = 1 in the expressions for the 
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FIG. 1. Photon loops with two external gravitational lines. 


/3-functions © Then these /3-functions © boil down 
to 


oIR 

J l,2,3 


(t - 1 ) = O, c). (18) 


10 (4tt) 


In order to find the anomaly-induced action one has to 
solve the equation 


= ~< T P = ^ + bE * + cUR ■ ( 19 ) 

V # OQ pis 

The solution of this equation has been originally found 
in ,23j and was discussed, e.g., in 01M1- Let us present 
only the final result for the covariant local form of the 
solution with two auxiliary fields [5H [35], 


r ind = Sclg^] + f d A xy/^g { * pA 4^ - * 


3c + 2b 2 
H 


1 


36 


■ V 


\f—b (£4 - - OR) 


S=b 


coC 2 


2 V-b 


</>C" 


}■ 


( 20 ) 


where 


A 4 = D 2 + 2 - -RD+ - R itl (21) 


is a covariant, self-adjoint, fourth-derivative, conformal 
operator [551 . 

The relevance of the auxiliary scalars ip and is based 
on the fact that the boundary conditions for these fields 
are equivalent to the boundary conditions for the two 
Green functions of the same operator A 4 in the non¬ 
local covariant form of Ti n d- The importance to have two 
auxiliary fields has been addressed in Refs. 
and [3], where the discussion of the role of conformal 
functional S c [g^, u \ can be also found. 

Let us make an important general observation. The 
dependence on the conformal factor in the Weyl-squared 
term is exactly the same as the one in the well-known 
logarithmic form factor for the massless quantum field^ 

T Ws = - l jJv^ 9 C apXT log (©) , (22) 


where T^s means T We yl-squared and /3i = (1607T 2 ) 2 for 
the photon. One has to remember that the term (22) 


comes from quantum corrections (loop of photon, in our 
case, see Fig. [T] for the illustration) independently on 
whether the higher derivative term is included or not 
into the classical vacuum action. Including this term is 
useful to make the procedure of renormalization more 
simple and regular. However, a finite term such as (22) 
will show up even if we set a 4 = 0. In what follows it 
will be shown that an arbitrary a\ will also come from 
fixing the initial condition for the soluti ons of the effective 
equations obtained by taking the term (201 into acco unt! 2 ! 
The same happens in the case of the quantum term ( |22| ), 
in the instant at which we fix the value of p. Therefore 
it is not critically relevant for our considerations that we 
started from the classical theory with higher derivative 
terms since the same terms emerge from quantum 
corrections in any case. This consideration shows also the 
deep relation between the anomaly-induced action and 
renormalization group. This relation will be extensively 
used in Sect. eh 

Starting from (201 one can derive the anomaly-induced 
effective equations for various cosmological epochs. Con¬ 
sider the effective action 


r — S vac + Ti n d j (23) 

where the classical action S vac was defined in (|T]) and the 
quantum correction T in d in (20). 

The equations for the auxiliary fields f/3 and <f> are 


yf~~9 

\f~~9 


A 40- 


V^b 


£4 - -OR - 


A 4 ^ - 


co 


.C 2 


3 

= 0. 


co 


2 


C 2 


2 

The conformal transformation of the metric 
9tiu(x) =gft V (x) e 2CT(x) , 


= 0 (24) 


(25) 


(26) 


gives the following transformation for the quantities of 
our interest (the reader can check [39] for details) 

V~9 A 4 = a/ 3 IA 4 , 

\f~~9 (^4 ~ ^ OR^j = y/—j} ^£4 — ^ □£ + 4A 4 cr^ . (27) 

For the sake of simplicity we consider only the spatially 
flat metric, 


g= V^a 2 {r]) = i) F e 2 " (,) 


(28) 


where 77 is conformal time, and use identities (27). Then 
Eqs. (24) and (251 become 

D 2 (f) = —2 'f—b0 2 <j , (29) 

□V = 0 , (30) 

where the d’Alembertian is the flat-space one. 


1 At the end of Sect. 6 we will discuss the IR limits of this corre¬ 
spondence in the asymptotically dS space. 


One see the details of this procedure in Eqs. fl44k, (|45|) and fl46|! 
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Solutions of Eqs. (29) and (30) can be found for vari¬ 


ous epochs of the cosmological evolution of our universe, 
i.e. separately during radiation, matter and during the 
de Sitter phase which is the asymptotic future for the 
ACDM universe. 


• During the radiation epoch, a(rj) oc rj and Eqs. (29) 
and (30) become 


4>W = 12 yf-brf*, 
^ = 0, 


(31) 

(32) 


where primes indicate derivatives with respect to confor¬ 
mal time. The solutions have the form 

ip = d 0 + di rj + d 2 r) 2 + d 3 r] 3 , (33) 

<t> = Co + Cl 1] + c 2 if + c 3 rf - 2 \T-b log (—') , (34) 

where do, ... d 3 , Cq, ... c 3 , are integration constants and 
Tjo is a generic reference time. 

• During the matter-dominated period, a(rj) oc rf and 
Eqs. (29) and (30) become 

<P (iv) =2AV^b V ~ 4 , (35) 

iP {iv) = 0 , (36) 

with solutions 

ip = do + do r] + d 2 rf + d 3 rf , (37) 

(p = c 0 + ci T] + c 2 rf + c 3 rf - 4log ( —) , (38) 

where do, ■■■do, Co, ■■■c 3 , are integration constants in 
general different from the corresponding ones in radia¬ 
tion. 

• During the de Sitter phase, we have a(ij) oc |^| _1 and 
Eqs. p9|) and ([30| become 


= -12 y/^bri 


-4 


(39) 

(40) 


ipW = 0 , 

with solutions 

ip = do + di rj + d 2 rf + d 3 r] 3 , (41) 

(p = c 0 + ci r/ + c 2 r/ 2 + c 3 rf + 2\f^b log ( — ) , (42) 

Here do, ...^3, Co, ...C3, are integration constants, 
in general different from the corresponding ones in 
radiation- and matter-dominated periods. 


IV. EFFECTIVE CORRECTIONS TO THE 
CLASSICAL ACTION 


• If a\ < 0 and k <C M p , the theory is stable under 
tensor perturbations, 

• If ai <0 and above some frequency k, which is 
comparable to M p , the theory starts to be unstable 
under tensor perturbations, 

• If ai > 0, the theory is unstable under tensor per¬ 
turbations V k. 

This result can be easily understood taking into ac¬ 
count the physical content of the higher derivative theory 
linearized around Minkowski space, as it was explained 
in Sect. [TTJ Our main interest is the dynamics of tensor 
modes in the late universe. According to [8], for ai < 0 
gravitational waves start an exponential growth only if 
their frequencies are close to the Planck order of magni¬ 
tude. For this reason in what follows we will restrict our 
analysis to frequencies satisfying k m 2 , this means 
that the Planck scale physics is beyond our considera¬ 
tion. One can assume that starting from the Planck scale 
the high energy gravitational theory passes through some 
qualitative change. For instance, in this regime the ap¬ 
propriate description may be (super)string theory, which 
is free of the problems of ghosts by construction (see sub¬ 
section 


VIB 5). 


In the late universe the cosmological background is 
varying quite slowly and one can define an effective co¬ 
efficient of the Weyl-squared term © , which takes into 
account semiclassical anomaly-induced effects ( |20| ), 

eff 


l {rf) = 01 + 


2 


[fPM ~ ^(v)] ■ (43) 


Using the results of Sec. m we obtain, in terms of phys¬ 
ical time t, 

radiation a = ap + A 3 1 1 ^ 2 + A 2 t + A 3 t 3 ^ 2 

+ i log (£)' (44) 

matter = ap + A 3 1 1 ^ 3 + A 2 t 2 ^ 3 + A 3 1 

+ T log G9- (45) 

de Sitter „ e // _ n c 


o, = ar 


■Ay e~ m + A 2 e~ 2Ht 


+ A 3 e 3Ht + ujHt , 


(46) 


In Refs. [TUI 


where to is a reference time, the coefficients A, are ar¬ 
bitrary integration constants and af contains both the 
classical ai and the constant part of the semiclassical 
contributions. In each epoch the A\ and of constant co¬ 
efficients may be different but for the sake of simplicity, 
we adopted the same notation for all the cases. 

The expressions (44), (45) and (46) come from the 


m and [8] tensor perturbations of the anomaly-induced action (|2()|). but they can be seen as a 


classical theory ( T]) around a FRW background have been 
investigated in details. The result is that the stability is 
fully determined by the sign of the coefficient ai of the 
Weyl-squared term and by the frequency of the pertur¬ 
bation k = I k I. It was found that 


local version of the renormalization group running for the 
parameter a 3 . Since the quantum electromagnetic field 
is strictly massless, one can trace this flow to the far IR. 
In the given theory the UV turns out to be more com¬ 
plicated, because starting from some high-energy scale 
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massive fields start to contribute to the running of a\. 
The contribution of massive fields breaks down the ele¬ 
gant form of the induced effective action, until the deep 
UV regime, when all the fields can be treated as massless 

again mm- 

Let us additionally comment on the expressions (44), 


(45) and (46). The complete solution for a\^ requires 
that the integration constants af and Ai ^,3 are deter¬ 
mined from experimental or observational data and the 
solutions for different phases are connected by requiring 
continuity. However, it is difficult to put this program 
into practise, because of the Planck suppression which 
makes impossible a direct observation of the effects of a\ 
or its running. The unique way to have an information 
about this coefficient is related to the evidence that the 
present universe is stable under tensor perturbations. We 
shall see in what follows that this is sufficient to obtain 
some physically interesting information. 

Our own existence shows that starting from the begin¬ 
ning of the radiation epoch and along the cosmological 
evolution up to present time, the sign of a\^ has been 
negative. If not, the tensor perturbations would have 
been of the tachyonic type with energy density of the 
Planck order of magnitude. Therefore, the proper fact of 
our existence puts strong restrictions on the values of of 
and j4i, 2,3 until the present time. 



t "tphM pi 


FIG. 2. Behavior of tensor perturbations in the case a\ < 
0, for different frequencies. The frequency k is measured in 
Planck units. 


and a massive spin-2 field h^J ' 1 , satisfying 

+ = 0. (50) 


V. TENSOR. PERTURBATIONS IN THE Since we are interested only in the tensor part of the 

ASYMPTOTICALLY DE SITTER PHASE perturbation h pv , let us consider 


According to the recent observations [321 S3] our un i" 
verse is now entering into a phase of its evolution domi¬ 
nated by the cosmological constant. The analysis of the 
time-dependent parameter a\^ (t) during the last phase 
of the evolution of the ACDM universe deserves a special 
consideration and will be explored in this section. 

Let us parametrize tensor perturbations around the 
de Sitter background g^ as 


9nv — 9^v T hf , 


(47) 


Since the variation of a\ in Eq. (46) is quite slow, it is 


possible to consider a constant a± in the equations for 
the tensor perturbations (see also discussion in [T5j and 

i)- 

The equations of motion in the transverse and traceless 
(TT) gauge can be cast into the form [8] [TO, 32] 


(□+y+^) (□ +^)V = 0, 


2A\ . 


(48) 


where □ is d’Alembertian operator on de Sitter back¬ 
ground and is the mass of the massive spin-2 mode 
m 2 = ~ (327rGai) _1 . Eq. (48) also agrees with 0H05]- 
The theory propagates a massless graviton hffi , satisfy¬ 
ing 



= 0 


(49) 


h^llS 


= h 


tens 

fllS 


( ° ) ' 


(51) 


where we used TT-gauge S^hf? = 0, d l hf 0 T = 0 and 
set a(t) = a(0) e Ht , with a(0) = 1. 

For the analysis of tensor perturbations it is more prac¬ 
tical to use physical time t. Furthermore, after we sep¬ 
arate the hJ 3 T - mode, there is no need to write indices 
and the field variable can be simply denoted as h = h(t). 
Then the Eq. (48) can be written as 


:: / 2 k 2 \ ■■ 

h + 6 H h + yllH 2 + m 2 ~\ - 2~J h 


+ ^6H i + 3H mi + 


2 Hk 2 \ ; /fc 4 m| k 2 




(52) 
h = 0, 


where we have used the background constraint A = 3 H 2 
and where dots indicate time derivatives. For a\ = 0, Eq. 

P reduces to the usual equation for General Relativity 
L) tensor perturbations around de Sitter background. 
The numerical investigation of Eq. ( [52] ) for both signs 
of ai has been performed in [8]. As we have already men¬ 
tioned, for frequencies k <C M p the case of a\ < 0 shows 
oscillating solutions without growing amplitudes. This 
behavior for three different frequencies is shown in Fig. 
[2] On the contrary, for a± > 0 one can observe a rapid 
growth of the metric perturbations for all frequencies, see 






















5000 


Hubble parameter in the pure de Sitter phase is 



t —^ph^pl 


FIG. 3. Dynamics of tensor perturbations for 01 > 0, for 
three different frequencies k = jkj. The plot for the lowest 
frequency has been rescaled to fit the figure. 


Fig|f] 

We want now to quantify which is the typical time 
scale at which anomaly-induced corrections to the clas¬ 
sical coefficient of the Weyl square term, ai, come into 
play. In pure de Sitter, we found that the behavior of 
is described by Ecp |46j), We focus our analysis on 
time-scales of order t ~ 1 [M p , since this is the expected 
scale of the instabilities. On this scales, the decaying 
exponentials in (46) vary slowly and can be treated as 
constant functions. Hence, we can parametrize Eq. (46) 
as 


al ff (t)=a 1 +u)H(t-t d s), (53) 

where tds is the time at which the density of matter 
and radiation become effectively negligible and a pure 
de Sitter phase starts. 

Since we are not assisting an overproduction of gravita¬ 
tional waves at present time, the nowadays value of a\ 
has to be negative (i.e., we are in the regime in which ten¬ 
sor perturbations are stable). As soon as * (t) crosses 
zero going to positive values, the gravitational waves en¬ 
ter the unstable phase (the massive ghost becomes tachy- 
onic). Then an exponential instability in the tensor sec¬ 
tor instantaneously shows-up. 

One can ask whether the stable phase, characterized 
by a negative will last for a long period of time. We 
can assume that the stable phase in which we arc living 
today will last at least until the beginning of the “pure” 
de Sitter phase of the expansion of our universe, hence 
we take a\ = a\^{tds) to be negative and of order one. 
In the ACDM model, with Ha = 0.7, the value of the 


3 Let us note that the numerical integration can not be performed 

for positive ai which are too much close to zero, since the growth 
of tensor perturbations is too fast in this case. 


H a = \ - ~ 0.8H 0 - 1.2 • 10~ 42 GeV. (54) 


Replacing H = H\ into (53) we find the time t q at which 


tachyonic modes emerge in the spectrum of gravitational 


waves, 


, _ l«il 
9 H A u 


~ 2.4 ■ 10 iJ yr = 2.4 ■ 10 3 4 bi. (55) 


Hence with the assumption |ai| ~ 1 the remaining time 
until the gravitational wave explosion is about thousand 
time longer than the time which already passed from the 
Big Bang. Of course, since we have no experimental data 
to define the present-day value of |5i|, the estimate given 
above has an ambiguity. According to Eq. (551 larger 


values of |ai| correspond to a longer stable phase, while 
smaller values of |ai| imply a shorter period of stability 
before the tachyonic explosion. 

Regardless of the quantitative side of this prediction, 
it looks interesting that our knowledge of quantum cor¬ 
rections to gravity is sufficient to know how the ACDM 
universe will end up. According to our analysis, there 
will be an instant gravitational waves explosion due to 
the tachyonic instabilities. 


VI. DISCUSSION OF THE ASSUMPTIONS 


We have started from a minimal set of hypotheses and 
arrived at the conclusion that the final destiny of the 
ACDM universe is not a peaceful infinitely long slightly 
accelerated expansion in the homogeneous and isotropic 
phase. On the contrary, the universe will end up in a 
strong tachyonic explosion of tensor perturbations. At 
first sight the intensity of this explosion (at least in the 
first order approximation) looks unrestricted, since at the 
instant when changes sign, the mass of the unsta¬ 
ble mode is infinite, according to Eq. (11). Then grav¬ 


itational waves of all frequencies will experience a fast 
exponential growth. However, since one can not trust 
the semiclassical approximation completely, the conser¬ 
vative estimate is that the tachyonic modes will have, at 
most, energy density of the Planck order of magnitude. 
Certainly, this is more than sufficient to provide a dra¬ 
matic effect on the geometry and matter contents of the 
universe. 

Another natural restriction comes from the fact that 
the instability which we have found corresponds to linear 
perturbations. It might happen that next orders in the 
perturbative expansion in will restore the stability. 
However, from the practical side both these restrictions 
do not matter too much, because even this “restricted” 
gravitational explosion should be capable to destroy the 
symmetry (homogeneity and isotropy) of the metric, and 
thus lead to the strong changes of the properties of space- 
time. 
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FIG. 4. One-loop diagram with electron loop and two external 
electromagnetic lines and the diagram representing the two- 
loop correction to the photon loop from Fig. [l] with the 
electron loop insertion. The diagram at the left is subject to 
decoupling in the IR. 


Since the result of our study looks so dramatic, it is 
interesting to consider the list of all the approximations 
which have been introduced in our analysis. Let us for¬ 
mulate various possibilities to avoid the instabilities in a 
the form of questions and answers. 


A. Completeness of the one-loop approximation 


Can we expect some qualitative changes in the result 
by taking higher-order loops into account? Formally, 
higher-loop contributions do not change the sign of the 
/3i-function j40| [41] , but this is not sufficient to draw the 
conclusion about the relevance of higher-loop terms. In 
the UV regime there will be higher logarithmic contribu¬ 
tions compared to the one-loop form factor (22) and this 


could produce strong effect on the running of a\. How¬ 
ever, the situation at low energies is quite different. Let 
us remember that second- and higher-loop corrections to 
the one-photon bubble include a loop of electrons or of 
other massive charged fermions, as shown in Fig. [4] Be¬ 
cause of the Appelquist and Carazzone decoupling theo¬ 
rem in QED, the contribution of such a loop is suppressed 
at least by a factor ( H 0 /m e ) 2 « 10^ 77 for the dynamics 
of the conformal factor. Since this is the part which gov¬ 
erns the running of aq , there are no chances that higher 
loops can change the result found for the conformal factor 
dynamics. 

At the same time, we have to remember that the goal 
of our study is not the dynamics of the conformal factor 
itself, but its interaction to gravitational waves. Indeed, 
the term in the effective action which we implicitly deal 
with is (22). At the two-loop order there will be non¬ 


local corrections to the photon propagator and photon- 
graviton vertex. Since we do not quantize gravity, the 
last type of diagrams can be ruled out. The simplest 
diagram at higher loops will be the one from Fig. fTlwith 
the polarization operator insertion shown in Fig. Hi due 
to the electron loop. 

Concerning this and further diagrams it is easy to see 
that the mass which defines the scale of decoupling will 
be the mass of the quantum fermion, while the typical 
energy of the external particle will be the one of the grav¬ 
itational wave. Finally, the relevant ratio which defines 
the effective cut-off for the higher-loop contributions to 


the running of a\ is (Sow Me) for a tensor mode with 
energy Sow- Of course this energy can be much big¬ 
ger than Hq. For instance, taking Sow — leV we have 
“only” ten-orders decoupling. Indeed, the gravitational 
wave production with corresponding energy density will 
have a very strong physical and geometrical effect. 

It is worth noting that we are dealing with the di¬ 
agrams which consist of the loops of matter fields and 
external lines of gravity. Therefore, the one-loop part is 
given by the contributions of free fields on a gravitational 
background. For this reason there is no dependence on 
the coupling constants at one loop order, and such de¬ 
pendence can be observed only starting from the second 
loop mm- Finally, the one-loop effects of photons are 
not suppressed by mass, since photon is massless. How¬ 
ever, all particles to which photon is coupled, are massive 
(and their masses are huge compared to the gravitational 
waves energy) and that is why the two-loop effects are in¬ 
deed strongly suppressed. 

The consideration presented above shows that the one- 
loop results ((441) , (45) and (46) are sufficiently robust 


in the framework of the semiclassical theory, at least for 
gravitational waves with frequency smaller than the elec¬ 
tron mass. The case of waves with higher frequencies de¬ 
serves an additional investigation, but it can not change 
the qualitative result which we obtained at the one-loop 
level. 


B. The possible role of quantum gravity 

What about quantum gravity (QG) effects, which have 
been neglected so far? It is not easy to give a definite 
answer due to the variety of existing models of QG. Let 
us consider a short list of the possibilities which have 
been better explored up to nowj^j 


1. Effective IR quantum gravity 


The standard effective framework for the IR effects of 
QG assumes that General Relativity (GR) is a universal 
theory of IR quantum gravity [49] (see also [50] for earlier 
work in this direction and m for an alternative opinion). 

The first thing to note is that GR is not a conformal 
theory, therefore a compact representation of quantum 
corrections in the form (20) is not possible. In this sit¬ 


uation one can make use of the renormalization group 
and of the perturbative /3-function /3\ for the parameter 
til. The jSi-function gives, in the case of massless theory 


4 Let us note that the problem of quantum stability of de Sitter 
space and the validity of semiclassical approximation has been 
extensively discussed in m and more recently in m and m- 
The difference between our conclusions and the ones of these 
works can be probably explained by the special role of tachyon 
tensor modes which we consider in the present work. 
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(e.g., quantum GR), the coefficient in the logarithmic 
form factor (22). Therefore, it can be still considered 
a good approximation in the present case. In this fre- 
mework the change of sign of the renormalization group 
equation (53) looks unclear, at least beyond the one-loop 
level. The reason is that higher-loop contributions of 
QG are beyond our control, because QG based on GR is 
non-renormalizable. Still we can draw some conclusions 
about the effective quantum gravity. At the one-loop 
level, the Weyl-squared quantum counterterm is known 
to be gauge-fixing dependent [52]. Then one can pro¬ 
vide any desirable value of 0 i by a special choice of the 
gauge-fixing parameters. This also means that the corre¬ 
sponding contribution vanishes on-shell and therefore can 
not be regarded as a physical effect. Furthermore, sim¬ 
ilar gauge-fixing dependence is expected for the leading 
logarithmic corrections at higher loops. Therefore, while 
the subject is not completely clear, one can suppose that 
the leading contributions of QG in this framework will be 
sub-logarithmic. As we have already seen in subsection 
VI A[ this means that the change of sign of 0i in © due 


to effective low-energy QG is very unlikely. 


2. Weyl conformal gravity 


In this case the form (20) of quantum corrections is 


available and the effect of QG is to increase the posi¬ 
tive value of /3\. The reason is that the contribution of 
conformal QG is positive [55H551 . exactly as the one of 
all matter fields |Tf| ). The expression for N v copies of 
massless vector fields has the form 




1 


(4tt ) 2 


'Nv 

v 10 


199 \ 
30 ) 


(56) 


where the last term is the QG contribution. It is easy to 
see that no cancellation is possible, and for N v = 1 the 
QG contribution shortens the period of time until the 
sign transition for of about one order of magnitude. 


3. General version of QG with fourth derivatives 


In this case, exactly like in the effective QG described 
in subsection |VI B T] the form of quantum corrections of 
the type ( 20 ), is not viable, because the original theory is 


not conformal, hence there is no anomaly to integrate. At 
the same time the /3-function for the parameter a i is well- 
defined, free of ambiguities [55] [53] [SB] and, according to 
well-verified calculations CHOCS EH] has the same positive 
sign as ©• The overall /3i-function for N v copies of 
massless vector fields has the form 


01 = 


1 


'N v l'66\ 

(4tt) 2 l To + "To"/' 


(57) 


In the physically interesting case we have one photon, 
N v = 1, and no cancelation can be expected. Moreover, 
in the far IR the ghost which provides some part of the 


QG contribution in (57), is supposed to decouple, and 
we come back to the situation described in subsection 

ivrBn 


4- Superrenormalizable QG 

The superrenormalizable models of QG have been first 
formulated in [58] on the basis of a polynomial action and 
shortly after that in [59] on the basis of an action of grav¬ 
ity non-polynomial in derivatives. In the polynomial case 
this model of QG has many (tensor and scalar) ghosts. 
In the non-local model of [59] there are no ghosts at the 
tree-level, however loop corrections lead to the emergence 
of an infinite amount of massive ghost-like states, corre¬ 
sponding to complex poles EDI- 

In both cases quantum corrections are well-defined, but 
the 0 i-function for the parameter ai can be modified by 
adjusting the terms of third- and fourth-orders in curva¬ 
ture tensor [60] - Therefore, it is not difficult to construct 
the theory with the cancelation of the photon contribu¬ 
tion in ©• However, this cancelation will take place in 
the UV region, where the quantum effects of ghosts are 
not suppressed. There is no systematic study of decou¬ 
pling of massive modes in higher derivative QG, but in 
principle one can expect that in the IR the massive ghost 
states will decouple and the low-energy situation will be 
covered by the effective QG, as described in subsection 
|VIB 1| After all, superrenormalizable models QG can 
not be expected to change the sign of 0i in the IR limit. 


5. String theory 

In (super)string theory the terms providing the /3- 
function for the parameter ai are usually removed by 
means of the Zwiebach transformation [331 eh m] ) - The 
procedure is ambiguous [63], but, by construction the 0- 
function for the parameter ai is zero. At the same time, 
string theory is not supposed to provide a significant cor¬ 
rection to the Quantum Field Theory results at low and 
very low energies, for otherwise we would observe such 
corrections in precision experiments, e.g., the ones that 
test QED and Standard Model calculations. Therefore 
using string theory to evaluate the /3i-function in the far 
IR is not reasonable from a conceptual point of view. 

In conclusion, we can see that QG can provide the 
change of sign of the /3-function for oi, but only at high 
energies and in the framework of superrenormalizable 
models of QG of [55j and [55]. On the other hand, the IR 
limit of the QG theory with a number of massive states 
(ghosts and normal particles) should be taken with great 
care HEEL] and it is expected that the sign of the 01 - 
function in IR would remain positive. Hence, we arrive at 
the conclusion that the main result concerning the pos¬ 
itiveness of the 0 -function for the parameter ai remains 
robust even if gravity is quantized. 
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C. Completeness of the anomaly-induced 
approximation 


As we have seen, there is practically no way to avoid 
the change of sign of ai within the anomaly-induced ef¬ 
fective action (20). However, we know that this effec¬ 


tive action is based on the artificial Minimal Substrac- 
tion (MS) scheme of renormalization. This scheme is 
always working perfectly well in the UV. On the other 
side, it is known that for the theory of massive quantum 
fields, this scheme is not working well in the IR, because 
of the decoupling theorem m, which was also obtained 
for semiclassical gravity in mm (see also J3] for the 
review and further references). In the flat space-time the 
masslessness guarantees that the main features of the UV 
will repeat in the IR, due to the UV-IR “duality” of the 
logarithmic form factor such as the one in (22). Is it true 


that the situation is the same in curved space-time? The 
question is not simple, especially in the asymptotically 
dS space, which has a natural IR cut-off scale H. 

For instance, the scalar curvature in the space with 
a(t) = Hot (with Hq = const) has the global scaling 
similar to the □ operator, namely R = — 6e~ 2<T ■ a" 2 
in the spatially flat case, but the derivatives here are 
with respect to conformal time rj, where dt = a(r])dr]. 
As a result R does not run with time, instead it is a 
constanl0 R = — 12H^. Similar consideration applies 
to the operator □ in the logarithm of the form factor in 
(22). For instance, let us consider □ acting on a scalar 


field <p. For the sake of simplicity we assume that this 
scalar depends only on time, Then a very 

simple calculation leads to the result 


□</? = (p + 3H 0 <p. 


(58) 


One has to note that in the last equation the depen¬ 
dence i p(t) corresponds to the fiducial (flat) metric and, 
therefore, time derivatives of ip are expected to be of the 
same order of magnitude in different epochs. Thus, Eq. 
(58) shows that the dependence of time drops out from 


the coefficients of the equation for perturbations. Sim¬ 
ilar calculation for the Weyl-squared term provides the 
result (here we kept the possibility to H be non-constant 
for the sake of generality) 

sj~ijC aPpT □ C af3pT = 

= ^C appT {df - 4 Hdt -2 H- 10 H 2 } C af)pT 
+ 8 ^§C afip 1 { Hd t - H + 2H 2 } C af)p *. (59) 

The dependence of space coordinates does not change 
the result qualitatively, if the wave vector is taken in the 
physical space, that is after the rescaling k -A k • e a . 

The expression (59) shows that in the asymptotic fu¬ 


ture, where the background becomes very close to de Sit¬ 
ter space, there will not be physical running of ai. The 


5 Authors are very grateful to A.A. Starobinsky for suggesting this 
example as an argument against IR running in de Sitter space. 


reason is that the extreme IR corresponds to dS and has 
its own scale H, such that the correspondence between 
logarithmic asymptotic ( |22| ) and anomaly-induced action 
(20) gets violated. Recently, explicit calculations of the 
curvature-curvature correlators on de Sitter background 
have been performed in j64j. It is remarkable that these 
calculations (performed by completely different method) 
have also shown the absence of leading-log corrections 
in far IR. Our interpretation of this result is a specific 
IR decoupling in the IR on de Sitter space background, 
where H plays the role of the natural minimal mass scale. 

Indeed, according to the same logic, the results for the 
radiation- and dust-dominated epochs, Eqs. (44) and 
(45) are valid and the same is true for the beginning of 
the transition period between matter and cosmological 
constant dominating epochs. However, since we are now 
in the situation when Hubble parameter is quite close to 
the future constant value H , this means the universe is 
entering the epoch of IR decoupling of all quantum effects 
in the higher derivative sector of the theory, including the 
ones of photons. Hence in the real situation the transition 
to the positive sign of a± is very improbable. 


VII. CONCLUSIONS 

The anomaly-induced effective action of vacuum is very 
efficient tool for a local implementation of the renormal¬ 
ization group running. According to this approach, the 
semiclassical effects drive the coefficient a± of the Weyl- 
squared term in the vacuum (gravitational) action to the 
positive side. It turns out that the one-loop contribu¬ 
tion is exact in this case, hence the future of the universe 
should be related to the tachyonic instabilities. Since the 
time-dependence of a\ is linear, one can easily calculate 
the time which remains until the change of its sign, that 
would be an instant when the massive ghost transforms 
into tachyonic ghost. 

The difference between the two types of ghosts is dra¬ 
matic. For the “usual” ghost there is a mass threshold 
and hence a possible situation when the ghost is harm¬ 
less. The situation can be understood in the spirit of 
considerations presented in [30]. Let us note that he ap¬ 
proach pursued in these papers by Simon and Parker and 
Simon is effectively equivalent to the one of [Sj. In both 
cases Planck frequencies are “forbidden” without clear 
explanation (we believe that one can find a solution of 
this problem, but this will be discussed elsewhere). Fur¬ 
thermore, in both cases a physically reasonable choice of 
initial conditions can be done, as a result the system does 
not fall into the run-away solutions. As a result higher 
derivative terms produce only tiny, Planck-suppressed 
corrections and can be simply ignored. For the gravi¬ 
tational waves this means that the dynamics of metric 
perturbations is essentially the same as in GR, as it was 
actually shown by Starobinsky in the second paper of [S]. 

As far as ghost becomes also tachyon, no mass thresh¬ 
old exist and no reasonable choice of initial conditions 
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can be done. Then an instant explosion of gravitational 
waves starts at all frequencies, independent on initial con¬ 
ditions. Our calculations show that in the framework of 
anomaly-induced effective action the transition to tachy- 
onic ghost is unavoidable in the far future of the ACDM 
universe. It is easy to show that the higher-loops correc¬ 
tions do not change this result, at least for the frequencies 
below the electron mass scale. 

However, in the dS-background case, which is a final 
stage of the ACDM universe, there is a qualitatively new 
type of IR decoupling which takes place even for massless 
fields. As a result, the anomaly-induced approximation, 
which is based on the MS scheme of renormalization be¬ 
comes non-appropriate in the IR, exactly as it happens in 
the massive theory. Therefore, there are no real chances 
of the tachyonic explosion of the universe even in the far 
future, since no change of sign of a i can be expected. 
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